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Abstract 

The problem of ^tt — > tttt is studied using the axial anomaly, elastic unitarity, analyt- 
icity and crossing symmetry. The solution of the integral equation for this amplitude is 
given by an iteration procedure. The final solution disagrees with vector meson dominance 
models. 
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One of the fundamental calculation in particle theory is the tt° — > 77 decay rate. It 
is a combination of Partial Conserved Axial Current (PCAC) and of the short distance 
behavior of Quantum Chromodynamics (QCD) [|IJ. The result of this calculation supports 
the concept of the color in QCD. This calculation is valid in the chiral world where the 
7T° is massless. Some correction has to be made, in principle, to take into account of the 
physical pion mass. It turns out that this axial anomaly formula is in very good agreement 
with the pion life time data, implying that the correction to the physical pion mass is 
very small. This result is expected because the pion mass is small. 

Another axial anomaly result is the process 771 — > or its analytical continuation 
7 — > 37r [0. This last process, even at moderate energy has a complex singularity and 
requires extra work and is postponed to a forthcoming study. 

The calculation of the process 77? — > nit is in itself interesting because not only it 
will be checked by future experiments which are being proposed at various accelerator 
facilities, but it will also be important for the calculation 7r° — > 77* H . 

Omitting the kinematical factor, the 7 — > 3tt and 77r — > tttt amplitudes are given as: 

F **(°) = 4^* = A C 1 ) 
where f n = 93 MeV and the zero in the argument of F 37T refers to the massless pions, and 
to the amplitude evaluated at the symmetry point s = t = u = 0. 

Using the KSRF relation || g p7T7T = m p /(\/2/ 7r ), it is straightforward to calculate the 
process 77? — > 7T7t |5|, [J using the vector meson dominance model ( VMD) : 

Fl MD (s, t, u) = \{^- + ^f- f + (2) 
2 rrip — s rrip — t m A p — u 

which in the chiral limit, s = t = u = 0, yields F^J 110 (0,0,0) = |A and is a factor 

3/2 larger than the anomaly Eq. ([!]). It is clear that the VMD model does not always 

work. In order to rectify this problem, a contact term is introduced to reduce the VMD 

contribution 0, |9|. The VMD equation can now be written as: 

m/»rn A rn? — ml m 2 n — ml m 2 n — ml 

Fg VMD (s,t,u) = ^ — ^ ^ + \ ; + -B- 1 - 1 3 

6 2 v m 2 p - s m 2 p -t m 2 -u 



and is seen to satisfy the anomaly Eq. ([!]) in the chiral limit. In Eq. (^) the smoothness 
assumption was made by setting F 37t (s — t — u — m\) = A. This assumption will be 
assumed to be correct in the remaining of this letter. 
Yet another solution was discussed ||: 

f mvmd (s t u) - x r 2 p- m " + m i~ m i , m 

^ ^ U >- 3 [ ml-s + m* p -t + m 2 p -u ) [ ) 

and is seen to satisfy the anomaly Eq. ([!]) but differs from the Vector Meson Dominance 
result by a factor of 3/2. 

In this article a different result on the amplitude based on the solution of an 
integral equation is given. The strong P-wave tttc scattering phase shifts is supposed to be 
known as given by the experimental data up to 1 GeV or higher, i.e. the existence of the 
p resonance at 0.77 GeV with a width of 0.151 GeV. We shall not make an assumption 
on the Vector Meson Dominance, but try to find a solution which is consistent with the 
constraint of the elastic unitarity, crossing symmetry and also of the low energy theorem 
Eq. (J]]). Our solution does not agree either with Eq. (Q) or with Eq. (|j). 

The process 77r — > tttt is a completely symmetric reaction, that is the same amplitude 
describes not only the reaction 77r° — > 7r + 7r _ , but also two other amplitudes involving the 
permutations of the pions. One assumes that the scattering amplitude can be represented 
by a single spectral function dispersion relation: 

F(s, t,u) = \+ [{a - ml) T -] + [s <- t] + [s - u] (5) 

J Ami (z — m^)(z — s — le) 

The spectral function can be evaluated, using the unitarity relation: n<j(s,t,u) = (n | 
73(0) I n)(n I 7T + 7r~) where the 73 is the isovector electromagnetic current operator and 
I n > is a complete set of the physical states; the summation over | n > is understood. 
Because of the G-parity, only the even number of pions intermediate states contribute. 
As one is interested in the low energy region (but including possible resonant effect) of 
the 77T° — > 7t + tt~ amplitudes where the inelastic effects are small, the summation over 
the intermediate states is truncated to the two pion (P-wave) states. The amplitudes 
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i^ I ,73(0) I 71+71 ) an d (7r + 7r I 7T + 7r ) are expanded into partial waves, and keeping only 
the lowest partial wave contributions (physically the P-state), one has: 

Tra(s) = F(s)e- i5{s) sin5(s) (6) 

where F(s) denotes the S-wave projection of F(s,t,u) (and physically corresponds to the 
P-wave when combining with the kinematical factor); 5 is the experimental P-wave irn 
phase shifts which can be obtained from the experimental data which shows that the 
phase shifts must pass through 90° at the p mass and its width is lhlMeV . 

Substituting Eq. (||) in the partial wave projection of Eq. (||) and keeping only the 
lowest partial wave, one has: 



(s-ml) f 00 F(z)e~ iS{z)) sm5(z) 
r (s) — AH / — - r\dz 



it J Am% ( z — m l)( z — s — ie) 



(7) 



The unitarity of the S-matrix requires that F(s) has the phase S |T(| and this should be 



exhibited in the solution for F(s). Because of the neglect of the contribution of higher 

partial waves, F(s) is not exactly equal to A at s = m%. Eq. (0) is a complicated integral 

equation; it is similar to, but more complicated than the Muskelishvili-Omnes (MO) type 

||, because the t and u channel contributions are also expressed in terms of the unknown 

function F(s). It should be noticed that the first term has a cut from 4m^. to 00 and 

the second one has a cut from to —00. This remark enables one to solve the integral 

equation by an iteration scheme as given below which converges very fast. 

The iterative and final solutions can be expressed in terms of the function D(s, 0), 

normalized to unity at s = and defined in terms of the phase shift 5 as: 

1 5 f°° 6{z)dz 

e M- / „ -77 — : — ( 8 ) 



D(s, 0) 7T J4mvr 2 z(z — S — it) 

Other functions D normalized to unity at s = Sq can be expressed in terms of the function 
D(s, 0) by the simple relation D(s,s ) = D(s, 0)/_D(s , 0). 



As remarked above, the Integral Equation (|7|) has both right and left cuts. This allows 
us to define an iteration procedure defined as follows: 



3 7r J4ml (z — ml)(z — s — ie) 

where is the value of the function F(s) calculated at the i th step in the iteration 
procedure; the Born term T % £ (s) calculated at the i th — 1 step is defined as 

Tt 1 \s) = -+- r F^\z)e- iS ^sm5(z){ 1 1 „ ln[l+ (a ~ 4m ^ ]-— J— } djg 

(10) 

where i > 1. The Born term is real for s > 0. 

The solution of the integral equation Eq. flj^) is of the MO type f8|: 

A , rrfi-i),^ , 1 s-ml D(z, ro%e**W sin 8(z)T%~ l \z)dz 

Uml (z — m%)(z — s - ie) 

(ii) 



3L>(s,m^) D{s,m*) it J a 



where it is assumed that the well-known polynomial ambiguity inherited in the MO inte- 
gral equation is absent ||. (The polynomial ambiguity could represent some incalculable 
inelastic effect occurring above the inelastic threshold). It is now straightforward to take 
out the imaginary part of the integral in Eq. (|TT|) and using e 4<5 sin<5 = p(s)N(s)/D(s,0) 
with ImD(s,0) = —p(s)N(s) to combine with the Tg -1 ^ Born term to show that FW(s) 
has indeed the phase 5 as required by the the unitarity of the S-matrix. 

One arbitrarily defines the convergence of the iteration scheme at the ith iteration step 
when | F« | / | F^ l ~^ | differs from 1 by less than 1% or so for the energy range from the 
2 7r threshold to lGeV. (Alternatively one can also consider the ratio | T$ \ / | Tg ^ |). 
Then combining the Tg ^ Born term in Eq. (|ll|) with higher uncorrected partial waves 
(for rescattering) from the t and u channels, one arrives at the final solution: 

F (i) (s, t,u) = £[{_J_(l + 3/^( S ))} + {(s - t)} + {(s <-> u)}] (12) 

where the function is the rescattering correction: 

/(*-!) ( s ) = g ~ m ' f°° D ^ Ky 5{z) stn5(z)T B - 1 (z)dz 
7i J Ami (z — ml)(z - s — ie) 



By projecting out the lowest partial wave from Eq. ( |l"2"D it can be shown that F^'(s) has 
the phase 5. 

In order to carry out the iteration procedure to find the solution of the integral equation 
one has to parameterize the function D(s,0) in terms of the experimental P-wave phase 
shift. Noticing that the phase of D(s, 0) is -5), one has [JTTJ, 0. 

_J_ = I (14) 

D(s,0) l-s/s R - ^± Tn {(s-Aml)H^{s)+2s/?>} V ! 

where H^(s) = (2 - 2^^Mn ^ ± ^^ L ) + m^J^f^ for s > 4ml; for other values 
of s, H n7T (s) can be obtained by analytic continuation. The partial wave amplitude is 
written as e i5 sm5/p(s) = N(s)/D(s, 0), where N(s) = (s - 4ml) / (96tt f%) . 

The p mass is defined as the vanishing of the real part of the denominator and is equal 
to sr in the narrow width approximation. Using y/s p =0.770GeV in Eq. flI3|), we have 
T p = 155MeV which is very close to the experimental value of 151.5 ± 1.5MeV. 

To make the iteration scheme coherent, the phase theorem for F^(s) has to be satisfied 
at every step of the iteration. One first guesses a solution for F^ which must satisfy 
this theorem; the following F^'(s) satisfies the theorem by construction. One first sets: 
F(°'(s) = (\/3)D^ 1 (s,ml); we then put this solution into Eq. ( |i~0D to evaluate the 0th 
order Born term: 

m n / n 2A . s n — ml . s — 4ml . , , 

3 s — 4ml s p 
To arrive at this relation one takes a 5-function approximation for the product 

e~ lS sin S(s)/D(s, 0) which is sufficiently accurate at this point. Eq. (|15|) can also be ob- 
tained from projecting out the lowest partial wave from the t and u channels contribution 
ofEq. §). 

Using this Born term, F^ is calculated, using Eq. ( PD with i — 1. 

One next calculates the next Born T B (s) term by using the solution F^(s) and then 
proceeds to calculate F^(s) etc. It is found that | F^(s) \ differs by less than 2.5% from 
| F^(s) | and | F™(s) \ differs by less than 0.4% from | F®(s) | in the energy range of 
2m,r to lGeV which shows indeed the very fast convergence of our iteration scheme. 
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If we were happy with an accuracy of about 1.5% in the amplitude, then the final 
solution given by Eq. flT2| ) could be simply expressed in terms of I^°'(s). The exact 
expression for J^(s) is a complicated linear combination of the Spence functions, but an 
approximate expression can simply be given as follows: 

/<0,(s) * L13 i4i% {Sp(1 + Ll 7r ) - Sp(1 - 3j f )] (16) 

where Sp(s) is the Spence function. With this approximation, the expression for (s, t, u) 
is accurate to within 2% compared to the final solution constructed from F^\s) in the 
energy range of 0.1 to 1 GeV. 

In order to compare these results with the previous pole models, one has to make a 
finite width correction for Eqs. This can be done by using the self-energy correction 

for the p propagator using the pmr coupling given by the KSRF relation. One then has 
to make the substitution (s p — m%)/(s p — s) by 1/D(s, m 2 ) and similarly for expressions 
with t and u variables. The lowest S- wave projection yields: 

and 

F 3 M J MD (s) = ^[—^ + n(s)} (18) 

where it is sufficiently accurate to use the pole terms for the contribution of the t and u 
channels. The phase theorem is violated in both Eqs. ( jT?| , |T8"| ) . 

Fig. 1 shows that the final solution lies between the two VMD models, the (CMVD) 
given by Eq. flTTD, and the (MVMD) model given by Eq. fllgp. It is useful to subject 
our result with future experimental test in the resonant region where the difference with 
other models are largest. 

Fig. 2 shows a similar plot but with only a low energy scale s < 0.3GeV 2 . The present 
experimental data at low energy Jl3|] confirms roughly the anomaly but is not sufficiently 
accurate to distinguish various models. With more accurate data, below s < 0.3GeV 2 , it 
will be experimentally difficult to distinguish our result with that of CVMD but it might 
be possible to distinguish our result with that of MVMD. For a more complete review of 
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the experimental situation as well as other theoretical models, two recent articles are to 
be consulted [TJ|, |T§ . 

It will be best to analyze future experimental data using Eq. (0), i.e. in terms of the 
energy and scattering angle instead of its S-wave projection as given by Figs. 1 and 2. 

In conclusion,the process 77r — > tttt is calculated using the low energy theorem, ana- 
lyticity, elastic unitarity and crossing symmetry which are fundamental conditions for a 
theory involving strong interaction. The final result shows that one effectively takes into 
account of the (unstable) p model in the s, t and u channels and their rescattering effect 
treated in a self-consistent way. It is important to put these results to experimental tests. 

Part of this work was done while the author was a visitor at the KEK National Labo- 
ratory and at the Hue University. The author would like to thank these two institutions 
for their hospitality and in particular, Prof. Y. Shimizu and his collaborators at KEK 
and Dr Nguyen trung Dan at the Hue University. 
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Figure Captions 

Fig.l:| F®(s) I 2 /A 2 clS cl function of the energy squared s(GeV 2 ). The solid curve is 
the solution of the integral equation after 3 iterations; the dotted curve is the solution 
after one iteration. The long dashed curve is from the modified VMD model (MVMD) 
Eq. (fL8l); the short dashed curve is from the VMD model with a contact term (CVMD) 
Eq. fln|). 

Fig. 2: | | 2 /A 2 as given by the Fig. 1 but with a smaller energy scale. 
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